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Abstract 

We analytically study stability of sliding lattice of Josephson vortices driven by a transport 
current in the stack direction in strong in-plane magnetic field. In contrast to recent findings we 
obtain that there are no diverse configurations of stable vortex lattices, and, hence, the stable 
sliding vortex lattice can not be selected by boundary conditions. We find that only the triangular 
(rhombic) lattice can be stable, its stability being limited by a critical velocity value. At higher 
velocities there are no simple stable lattices with single flux line per unit cell. Oblique sliding 
lattices are found to be never stable. Instability of such lattices is revealed beyond the linear 
approximation in perturbations of the lattice. 

PACS numbers: 74.25.Fy, 74.80.Dm 
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I. INTRODUCTION 



Strong anisotropy of layered superconductors, such as high-T c superconductors, results in 
a number of specific features. In particular, magnetic field parallel to the conducting layers 
induces coreless vortices similar to fluxons in conventional Josephson tunnel junctions 1 ' 2 . 
In contrast to the fluxons in conventional Josephson junctions, the current and magnetic 
field of a Josephson vortex in layered superconductors are spread over many conducting 
layers. The transport current flowing in the stack direction forces the vortices to slide along 
superconducting layers. Since in contrast to standard Abrikosov vortices the Josephson 
vortices are coreless, their motion does not involve perturbations of the amplitude of the 
superconducting order parameter and corresponds at low temperatures to an underdamped 
regime. So the velocity of such vortices, in principle, can be quite large. Many dynamical 
properties of Josephson vortices in layered superconducting crystals are similar to those 
of fluxons in artificially prepared stacked Josephson junctions (see, e. g., Ref. 3). The 
Josephson vortices can be arranged in a lattice, however, possibility of multiple metastable 
states corresponding to vortex configurations which do not form any lattice is predicted as 
well at small magnetic field 4 . At larger magnetic field the Josephson vortices form a lattice 
and, hence, in the flux-flow regime the coherent motion of the vortex lattice is expected. 
Such a regime must induce coherent electromagnetic radiation from the uniformly sliding 
lattice. The flux-flow regime of this lattice was observed in BSCCO mesa structures 5 ' 6 ' 7,8 ' 9 . 
However, the regime of coherent sliding of the lattice was found to be limited by a maximum 
voltage and, hence, by a maximum lattice velocity 6 ' 7 ' 8,9 . Above this voltage either a different 
regime was observed 7 ' 8 or no stable I-V curves were found in a close range of currents and 
voltages 6 . In addition, above the maximum velocity a broad-band non- Josephson emission 
in the microwave region was observed 7 ' 8 . 

The upper limit for the velocity of the coherent flux-flow regime was obtained theoretically 
in our studies of stability of the vortex lattice 10 ' 11 , and was related to the interaction of the 
vortices with Josephson plasma modes excited by moving vortices. A regular motion of 
Josephson vortices in a form of the rhombic lattice, usually called the triangular lattice, 
was found to be stable at velocities up to a critical value v cr ~ cs / (2X\\y/e) , where s is the 
period of the crystal lattice in the stack direction, An is the London penetration depth, e is 
the dielectric constant in the stack direction. This critical velocity was identified with the 
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experimentally observed velocity limiting two regimes of the vortex motion 6 ' 7 ' 8,9 . We have 
found that at higher velocities sliding of a regular vortex lattice is frustrated due to a growth 
of fluctuations of the vortex lattice. The instability is induced by interaction of the vortices 
with the Josephson plasma mode. The similar results were obtained both in the limit of 
large magnetic fields, when the distance between the vortices is smaller than the size of the 
non-linear region near the vortex center, and for the case of lower magnetic fields 10 ' 11 when 
the vortices are separated by distances much larger than the size of the non-linear region. 

However, qualitatively different conclusions on the stability of moving vortex lattice were 
made in recent papers by Koshelev and Aranson 12 ' 13 . According to their calculations a 
stable vortex lattice motion can occur not only in a form of the triangular lattice, but 
also in a form of various oblique lattices (see fig. 1), a particular experimentally realized 
lattice being uniquely selected by boundary conditions. Regimes of stable motion of various 
oblique vortex lattices were found also for velocities exceeding the critical value v cr . The 
existence of a set of various stable vortex lattices for a given sliding velocity would lead 
to interesting physical consequences. Since a particular vortex lattice is uniquely selected 
by boundary conditions, the radiation spectrum must also strongly depend on boundary 
conditions. Specifically, in contrast to other vortex lattices, the triangular one does not emit 
at the main radiation frequency related to the time needed to shift the lattice by its period. 
This happens since the emissions by adjacent vortex rows, which in the triangular lattice are 
shifted by the half-period, compensate each other. Therefore, in case of the triangular lattice 
the radiation is generated at even harmonics only. The oblique lattices, in which the shifts 
between the adjacent vortex rows are not equal to the half-period, can radiate both at even 
and odd harmonics as well. Furthermore, the intensity of the radiation maximum predicted 
at the Josephson plasma frequency 10,11 would be much larger if there exists a stable vortex 
lattice arrangement at v > v cr . 

To eliminate the discrepancy and to resolve the problem we reconsider stability of sliding 
vortex lattices and find that the oblique vortex lattices are not stable. Instability of oblique 
lattices is found out from the analysis of equations of motion for perturbations of the lattice 
beyond the linear approximation in the perturbations (Sec. IV). We illustrate this statement 
by a particular case of the equilibrium state, i. e., of the lattice at rest (Sec. III). In the 
equilibrium state the results are especially transparent, because the stable vortex arrange- 
ment can be selected as that corresponding to a minimal energy. We conclude that oblique 
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lattices correspond to a saddle point of the energy functional with respect to perturbations 
of a uniform lattice, and only the triangular lattice corresponds to a minimum. This is 
revealed when the terms beyond the quadratic ones are kept in the corresponding energy 
expansion. Therefore, to study stability of sliding vortex lattices by means of equations of 
motion one should proceed beyond the linear approximation with respect to perturbations. 



II. MAIN EQUATIONS 

Electro dynamic properties of layered superconductors can be described in terms of gauge 
invariant potentials which can be treated as superconducting momentum in layer n, p n , and 
gauge invariant phase difference between layers n + 1 and n of the superconducting order 
parameter, (p n . In general case electric and magnetic fields in the superconductor depend 
also on gauge-invariant scalar potential \i n = {l/2)d t Xn + &n, &n is the electric potential in 
the n-th layer. We do not take into account here the components of the electric field related 
to /i n since we consider low temperatures when the effects of branch imbalance related to 
H n can be ignored. 

Equations for the moving vortex lattice are derived similar to our previous paper 14 , 
namely, substituting the expression for the current density j presented as a sum of su- 
perconducting and quasiparticle currents to the Maxwell equation 

4tt 1 
V x H = — j + -dtD. 

c c 



Then we find equations for p n and ip n which have a form: 



d x p n+ i - d x p n 4vr 



j c sin ip n + —d t ip n 
2s 



+ 7T~ d tt<Pn 



2cs 



2s 2 



(d x <Pr. 



C 4:71 

Oxipn-i) + — (p n+1 + p n -i - 2p n ) = 



-p n + G\\dtPr, 



(1) 

(2) 



Here Ay and Aj_ are screening length for screening currents flowing in the plane and stack 
directions, respectively, cry and <t_l are quasiparticle conductivities for related directions. 

Excluding superconducting momenta from Eqs. (1-2) one readily obtains the equation of 
motion for superconducting phases ip n 



(dlip n + uj r d t ip n + u 2 p sin y? n ) = ^ ^1 + ^d^j d 2 xx ip n . (3) 
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Here d 2 f n = fn+i + /n-i — 2/ n corresponds to the discrete version of the second derivative. 
Q r = 4man, u r = A7ra±/e are relaxation frequencies. Q p and uj p are plasma frequencies 
for the directions along and perpendicular to the layers, Q p = c/Aii is much larger than all 
typical frequencies of the considered problem, u p = cj y/e\± <C fl p . 

Studies of the vortex dynamics both in our papers and in papers by Koshelev and Aranson 
are based upon such equation. 

We limit our study by the case of large magnetic field considered in Refs. 12,13 . In high 
magnetic field the cores of Josephson vortices, i. e., non-linear regions at the centers of 
Josephson vortices, strongly overlap and solution of Eq. (3) can be found perturbatively as 
a sum of a linear term and a small oscillating correction: 



(0) 



Y n + i) n {Y n ). 



(4) 



Here Y n = q x — oj^t + ojq = qov, qo = 2tt/S\\, S\\ is the period of the chain of vortices 
along the superconducting layers, and x specifies the type of the vortex lattice (see fig. 1). 
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FIG. 1: The configuration of the vortex lattice is defined by % = 27r<5*S r n /'S'li • The rectangular 
lattice corresponds to x = 0, the oblique lattices are defined by < x < and the triangular 
lattice corresponds to x = 7r - 



, iuj^K 2 (ujn, x) 

Here K(oo,q ± ) 2 = 1 - iVt r uj/uj 2 p + 4(A,|/s) 2 sin 2 q ± /2, c 2 (uj) = c 2 (l - z^cu/cj 2 )/^. The 
denominator of G(uiq, qo, x) does not become zero since u = q v and vortex velocity is small 
enough if below the critical value v cr mentioned above and discussed below in Sec. IV. 
Hence, for stable vortex lattice |-0 n (F)| <ti 1. 
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III. STATIONARY VORTEX LATTICE 



At first we consider the vortex lattice at rest as the most simple and physically transparent 
case. Stability of a stationary vortex lattice is determined by its free energy. Minimum of 
the free energy as function of x corresponds to the stable vortex lattice. The energy of 
the vortex lattice can be considered as a sum of the energy of magnetic field induced by 
superconducting electrons, 

EW = J2 I dxH 2 /8n, 
of the kinetic energy of superconducting currents 

and the Josephson energy of n Junctions. 

n ^ 

Then the total energy can be readily presented as 



E (tot) 



dx 



^ i^f - <p'n(Pn + l - Pn) + ( P n + 1 ~ Pnf + ? f ) + |(1 " COS Vn ) 



(6) 

The variation of the energy functional (6) over p n and ip n yields Eqs. (1-2) with dissipative 
terms equal to zero. 

To determine the stable vortex configuration we substitute expressions for p n (Y n ) and 
<^4°) (4-5) to equation for the energy density (6). After some algebra the expression for the 
energy density of the vortex lattice reads 

{tot) _ J_ / 2ZW(f ) 2 (sin(f ) - I) 2 _ / 2 sin 2 (f) 1 1 1 \ 

~8ns{ ~ XWo Alg 2 + 4 % 8\Wo ^l) ' {) 

From Eq. (7) one can see that the energy minimum corresponds to x — n - It means that 
in accordance with previous results (see, e. g., Ref. 2), only the triangular lattice is stable 
(fig. 2). The difference between the energies of the triangular and the rectangular lattices is 
small and proportional to q^ 2 . 



IV. STABILITY OF SLIDING JOSEPHSON VORTEX LATTICE 

Now we consider a sliding vortex lattice. We adopt here the method used in Ref. 15 
to study a system of two coupled Josephson junctions. The solution of (3) for <p n can be 
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presented in a form 

<p n (x,t) = <p<®(Y) + 6<p n (x,t), (8) 

where ^\Y) given by Eqs. (4-5) describes the uniform motion of the vortex lattice with 
the velocity v, and 8<p n (x, t) corresponds small perturbations of the uniformly sliding lattice, 
\6<p n (x,t)\ < 1. 

The vortex lattice is stable if all possible perturbations Sip n (x,t) decrease with time and 
there are no increasing solutions for 8ip n (x,t). If in the linear approximation in 8(p n (x,t) 
there are solution which neither decay nor increase, then the linear approximation in per- 
turbation is not sufficient and one should take into account higher order perturbations in 
6ip n (x,t). 

The equation for 5(p in the linear approximation we derive substituting Eq. (5) into Eq. (3) 
and performing Fourier transformation over x, t and n. It reads 

F(uj, q h q ± )5tp(u, q\\, q±) + f(u + u , q\\ + q Q , q± + + <7|| + 9o, q± + x) + 

f*(u - u , q\\ - g , q± - X)8<p(w - u , q\\ - q , q±~x) = °- ( 9 ) 

Here q\\ and q± {\q±\ < 7r) are the wave vectors in the directions parallel and perpendicular 
to the superconducting layers, respectively, F(u, q\\, q±) = c(u) 2 qj — K 2 (u, q±)(uj 2 + iuui r — 
iu 2 ^g(uj ,q ,x)), f(u,q\\,q±) = K 2 (uj } q ± )u;l/2. 

The values of the perturbation with shifted arguments, 5ip{u) ± ojq, q\\ ± q , gj_ ± x), are 
small compared to 6ip(u),q\\,q±) with small values of the arguments. So using perturbation 
approach with respect to S\\/\j and neglecting 5ip(u±2uJo, q\\ ±2q , qj_±2x), we express the 
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shifted values in terms of 5<p(u,q\\,q±) and find a simple equation of motion for 5ip which 
determines the eigenmodes of the vortex lattice. 

fjpf \ f*(w,q\\,q±)f(w + wo,q\\ + qo,q± + x) 

* [u, q\\, q±) — ■ ■ ■ — r 

V F(uj + u , q\\ + g , q± + X) 

/(w, gy, g±)f*(u - uo, q\\ - go, g± - x) 



8<p(w,q\\,q±) = Q- ( 10 ) 



F(u - u , gy - g , q± - x) 
In this expression two last terms originate from interaction of the first harmonic of oscillating 
field induced by the sliding lattice with the Josephson plasma mode. 

As it follows from Eq.(5) Q(u ,q ,x) oc (S\i/\j) 2 is small, and one may neglect terms 
with G(uJo,qo,x) in the expressions for F in the denominators of the last terms of Eq.(lO). 
Such simplifications allow us to calculate from Eq.(lO) the dispersion relation defining the 
spectrum of collective oscillations 

2 • C 2 (cj) 2 
U = -lUJUJ r + — r + U) 

K 2 (uJ,q ± ) 

x 3 (iG(u Q , g , x) - ^G(w + w , q\\ + go, q± + x)~ ~ ^o, q\\ ~ go, q± ~ x)) • (H) 

Consider first the triangular vortex lattice, \ = n - 111 the long-wavelength limit, for 
the lattice in rest we obtain the sound-like spectrum with the velocities cy = c/^Je and 
c_i_ = (AS , j|/v / 8e7rA^)c in the in-plane and in stack directions, respectively. These velocities 
coincide with ones obtained by Volkov 16 . For the sliding triangular lattice the imaginary 
part of the frequency of the eigenmodes, which determines the decrement of damping, is 
positive for small velocities v < v cr ~ (s/2 v / eA||)c. Thus at such velocities the triangular 
lattice is stable. 

As it follows from Eq. (11) at larger velocities for any value of x there are values of gx 
for which the imaginary part of the frequency of the oscillations becomes positive. This 
means that fluctuations increase with time resulting in an instability of the lattice at such 
velocities. The value of v cr is about 3 x 10 5 m/c for typical parameters Ay/s = 100 and 
e ps 20. This estimate of the critical velocity coincides with the velocity limiting flux-flow 
regime in Refs. 6,7. 

Now we consider in more details the case of arbitrary x an d v < v cr . The sign of 
the imaginary part of u>, which determines the decrement of damping, sets conditions for 
stability of the vortex lattice. A positive damping factor for a given x an d for all values 
of q\\ and g_|_ means the stability of the vortex lattice. While a negative damping factor, 



in other words, a positive increment of growth, manifests the instability. Analytical and 
numerical calculations show that the damping factor is not negative for the whole region 
of values of x f° r a given value of the vortex velocity including zero velocity. This result 
was interpreted in Refs. 12,13 as an evidence of the stability of these lattices. However, 
the existence of a set of different stable lattices at zero velocity apparently contradicts to 
the studies of the equilibrium lattice 2 and to the results of the preceding section according 
to which only the triangular lattice should be stable. The point is that though for the 
oblique lattices the damping factor is positive for finite values of q, it is equal to zero 
for q = 0. This indicates that the linear approximation on Sip does not give an ultimate 
answer on the stability problem, and further study is needed. The situation can be clarified 
considering the particular example of the equilibrium state when the stable solution can be 
selected using the energy considerations. If the series expansion of the energy with respect 
to perturbations near the extremum does not contain the quadratic term, then contribution 
of the potential energy to the linearized equation of motion becomes zero suggesting a zero 
damping coefficient in a corresponding equation of motion. Therefore in order to study 
stability in such a case one must take into account the next order perturbations in the 
energy expansion in perturbations. If the third order term in perturbations in the expansion 
of the potential energy is non-zero then the extremum of the energy corresponds to a saddle 
point and, hence, to an instability. 

These considerations cannot be directly applied to the nonequilibrium case the stability 
of which should be studied by means of equations of motion, but they point out that to find 
whether solutions with zero damping factor in the linear approximation are really stable, one 
must include in the equation of motion for the phase the next order terms in perturbations 
5<f. Keeping in Eq.(2) linear and quadratic terms in perturbations we find the equation for 
Sip which has the form 



c(uj) 2 q 2 - K 2 (uj, q±)(uj 2 + iujuj r - iui 2 p $ (p(uj, q\\,q±j)) <fy(o) + 
K (+i) Sl P(+i) + ^( 2 _i)<fy(-i)] + ^ 9? (p(w, q\\,q±j) K 2 (uj, q±) [Sip 2 



-o) 



V 

2i 



K 2 

A (+i) 



K 2 



= 0. 



(12) 



Here 5ip (0) = 5ip(uj, q\\,q±), 5</?(±i) = 5ip(uj±u , q\\±q , q±±x), Kf±i) = K 2 (u±ujq, q±±x), and 
[<V 2 ±1 ) denote the Fourier transforms of [5ip n (t, x)] 2 with arguments (lu±lu , q\\±q , q±±x)- 
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Eq. (12) can be studied by perturbation approach in Sip and (S\\/\j) 2 similar to previous 
study of Eq. (9). Furthermore, we limit our analysis of the nonlinear case by the study 
of stability with respect to the uniform perturbations, since the stability of the state with 
q = only was not ascertain in the linear approximation in 5ip. As the first approximation 
for Sip we choose a solution corresponding to the perturbation of the parameter \ describing 
the tilt of the lattice: 

V 1J = \m5x + 5(0(^0, qo,xW Y ) = ^^((l + ^(x)sin X )e^). 

Here Gi(x) — dG(^o, Qo, x)/<9( sm x/2)- Terms with e 2lY and e~ 2lY give small corrections to 
(5(p^) 2 and should be neglected. Finally, the equation for vortex lattice perturbation (12) 
in the second order can be rewritten as 

(uj 2 + iuu r )6<pM = 2 3fc&( X ) sm(x)(5 X ) 2 . 

Since the right-hand size of this equation does not depend on u and ^tQ\{x) 7^ for any 
value of X) there are always increasing perturbations 5(p^ unless the factor before (Sx) 2 
is equal to zero, i. e., when x is a multiple of n. As it was found above in the linear 
approximation in lattice perturbations the rectangular lattice, x — 0, is not stable due to 
perturbations with finite q±, so there is only one simple vortex lattice with one flux quantum 
per cell which could be stable, namely, the triangular one, x — 71 ■ 

V. CONCLUSIONS 

In the frame of a simple theoretical model we studied stability of sliding lattices of 
Josephson vortices induced by strong magnetic field applied parallel to conducting layers of 
layered superconductor and driven by a transport current in the stack direction. We found 
that in contrast to recent theoretical predictions 12 ' 13 the type of the stable sliding lattice 
is not selected by boundary conditions, but there is only one possible stable vortex lattice 
arrangement of sliding Josephson vortices containing single flux line per unit cell, namely, 
the triangular lattice. In other words, an experimentalist cannot have an influence on the 
structure of the sliding vortex lattice manipulating with experimental conditions. Sliding 
of vortices in the form of the regular lattice is found to be stable up to the critical velocity 
only. At higher velocities there is no sliding regime with the vortices coherently arranged in 
an oblique, triangular or rectangular vortex lattices. 
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